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PDFs for the LHC 

Amanda Cooper-Sarkar 

University of Oxford, Denys Wilkinson Bdg, Kcblc Rd, Oxford 0X1 3RH, GB 

The PDF4LHC group has benchmarked six modern PDF sets and used them to make 
predictions for W,Z prdocution and Higgs production at the LHC. The reasons why pre- 
dictions differ are examined and recent updates to the PDFs and their predictions are 
presented. 



1 Introduction 

PhI The PDF4LHC group have benchmarked modern PDF sets from the groups; ABKM, CTEQ, 

p ^i GJR, HERAPDF, NNPDF, MSTW in terms of the predictions for basic LHC cross-sections [T|. 

D ' Fig [T] shows the comparison of predictions for the W'^ cross section at 7 TcV. The PDF4LHC 

r^ '. recommendation [2 is to use the envelope of the MSTW2008, CTEQ6.6 and NNPDF2.0 pre- 

dictions and add the uncertainty due to a^i^Mz) in quadrature. However, as can be seen from 
the figure, this may not be sufficient to cover PDF uncertainties. 

^ [ There are several reasons why the PDF predictions differ: firstly they arc based on different 

data sets. For example, all the 2010 PDF analyses illustrated, bar the HERAPDFl.O and the 
NNPDF2.0, do not use the recently combined inclusive cross section data from HERA-I [3] 
l/-v . which are up to three times more accurate than the separate HI and ZEUS data sets used 

by previous PDF analsyses. These combined HERA data are also shifted in normalisation 
by ^ 3% with respect to the previous HERA data, and this partly explains the higher cross 
section of the HERAPDF wrt CTEQ and MSTW. Secondly, PDFs use different central values 
of as{Mz), the effect of this is illusrated on the figure. Some groups (HERAPDF, CTEQ, 
NNPDF) adopt a central value of asiMz) inspired by the PDG value and others (ABKM, 
GJR, MSTW) fit as{Mz) simultaneously with the PDF parameters and use their best fit value. 
Thirdly, the PDF analyses differ in the schemes used to account for heavy quark production. 
For example, the NNPDF2.0 used a zero-mass variable-flavour number scheme (ZMVFN) and 
this explains why the NNPDF2.0 predictions lie lower than CTEQ, MSTW, HERAPDF all of 
which use general mass variable flavour number schemes (GMVFNs). The value of the charm 
and beauty masses also differ between the PDF analyses. Lastly, PDFs differ regarding choices 
of PDF parametrisation and theoretical/model prejudices which are imposed. Whereas this 
latter source of difference represents legitimate, irreducible differences in approach it may be 
possible to achieve greater concordance on the flrst three reasons for PDF differences. 

2 Heavy Quark Schemes 

Let us flrst consider heavy quark production. The ABKM and GJR groups use Fixed Flavour 
Number (FFN) treatments, HERAPDF, CTEQ and MSTW use GMVFN and NNPDF2.0 used 
ZMVFN. However even within GMVFN schemes there is not complete agreement. Predictions 
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Figure 1: Predictions for the W+ cross section at the LHC at 7 TeV from various modern 
PDF sets as a function of as{Mz)- The cross section for each PDF set is shown at the value 
of as{Mz) used by that group. The vertical error bars represent the 68% uncertainty on the 
predictions, the horizontal error bars represent the 68% uncertainty on as{Mz) considered by 
each group. The right hand plot shows the predictions for PDF sets available in April 2010, 
the left hand side shows the updates from CTIO and NNPDF2.1 available for April 2011. Plots 



from G.Watt http://projects.hepforge.org/mstwpdf/pdf41hc/2010/ 



for F2 differ between schemes [4] and the choice of scale within a scheme affects predictions. 
The value of the charm mass can also affect predictions, HERAPDF, NNPDF and MSTW now 
provide PDFs at different charm mass values so that the effect of this can be evaluated. HI 
and ZEUS have recently combined their data on ^2'^ [5], and these data can help to reduce 
the uncertainty on PDFs coming from the choice of scheme and the value of the charm mass. 
These data have been input to the HERAPDF fit together with the inclusive data which were 
used for HERAPDFl.O. The x^ oi this fit is sensitive to the value of the charm quark mass. 
Fig. [2] compares the x^j a-s a function of this mass, for a fit which includes these data (top 
right) to that for the HERAPDFl.O fit (top left). However, it would be premature to conclude 
that the data can be used to determine the charm pole-mass. The HERAPDF formalism 
uses the Thorne- Roberts (RT) variable- fiavour-number (VFN) scheme for heavy quarks. This 
scheme is not unique, specific choices are made for threshold behaviour. In Fig. [2] (bottom left) 
the x^ profiles for the standard and the optimized versions (optimized for smooth threshold 
behaviour [6]) of this scheme are compared. The same figure also compares the alternative 
ACOT VFN schemes and the Zero-Mass VFN scheme. Each of theese schemes favours a 
different value for the charm quark mass, and the fit to the data is equally good for all the 
heavy-quark-mass schemes excepting the zero-mass scheme. Each of these schemes can also 
be used to predict W and Z production for the LHC and their predictions for W~^ are shown 
in Fig. [2] as a function of the charm quark mass (bottom right). If a particular value of the 
charm mass is chosen then the spread of predictions is as large as ~ 7%. However this spread 
is considerably reduced ~ 1% if each heavy quark scheme is used at its own favoured value of 
the charm quark-mass. Furher details of this study are given in ref. [?]• 

A further relevant question is what is the charm mass which is being used in the calculations? 
Is it the pole mass or the running mass? The running mass is independently measured and the 
ABM group have considered using this mass in their FFN calculations [5] • They obtain a value 
of 77T,c(mc) = 1.18 ±0.06 GeV in agreement with the PDG value. This should be improved after 
the input of the combined HERA charm data. 
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Figure 2: The x^ of the HERAPDF fit as a function of the charm mass parameter m'^"'^^^ . 
Top left; using the RT-standard heavy-quark- mass scheme, when only inclusive DIS data are 
included in the fit. Top right; using the RT-standard heavy-quark-mass scheme, when the data 
for ¥2'^ are also included in the fit. Bottom left; using various heavy-quark-mass schemes, when 
the data for ¥2'^ are also included in the fit. Bottom right: predictions for the W^ cross-sections 
at the LHC, as a function of the charm mass parameter to™°'''^', for various heavy-quark- mass 
schemes. 



3 2011 updates 

There have been recent updates to the PDFs which address some of the issues we have raised. 
NNPDF2.1 [5 updates NNPDF2.0 to use a GMVFN. CTIO |2] is an update of CTEQ6.6 
which includes the combined HERA data in addition to giving special consideration to the 
Tcvatron W asymmetry data (CTIOW). ABMll [10] is an update of ABKM which includes the 
combined HERA data. HERAPDFl.O has been updated to HERAPDF1.5 [U] by including 
further preliminary combined data from HERA-H running. Preliminary HI data on NC and 
CC e/^p and e~p inclusive cross-sections and published ZEUS data on NC and CC e^p and 
CC e+p data, from HERA-H running, have been combined with the HERA-I data to yield 
an inclusive data set wih improved accuracy at high Q^ and high x |12| . This new data set 
is used as the sole input to a PDF fit called HERAPDF 1.5 which uses the same formalism 
and assumptions as the HERAPDFl.O fit. Fig. [3] (left) shows the combined data for NC e^p 
cross-sections with the HERAPDF1.5 fit superimposed. The parton distribution functions from 
HERAPDFl.O and HERAPDF1.5 arc compared in Fig. [3] (right). The improvement in precision 
at high X is clearly visible. 

Fig. [1] right hand side, shows the predictions lot the W+ cross section at the LHC at 7 TeV 
updated to show the predictions of the new PDFs from CTIO and NNPDF2.1. The CTIO and 
CTEQ6.6 predictions are very similar and the HERAPDF1.5 (not shown) and HERAPDFl.O 
predictions arc very similar. The NNPDF2.1 prediction has moved up significantly beacusc of 
the use of the GMVFN scheme. The early LHC data agree well with all the predictions. 

It is also interesting to look at the predictions for the ratios of cross sections Z/{W^ + 
W-) and W-/W+. Fig. H shows 2010 and updated 2011 predictions. The Z/W ratio is 
very consistently predicted because fiavour dependence almost cancels in the ratio, but the 
W~^/W~ ratio predictions are quite discrepant - this measures a difference in the u and d 
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Figure 3: Left: HERA combined data points for the NC e^p cross-sections as a function of 
Q'^ in bins of x, for data from the HERA-I and II run periods. The IIERAPDF1.5 fit to these 
data is also shown on the plot. Right: Parton distribution functions from HERAPDFl.O and 
HERAPDF1.5; xuy, xdy^xS = 2x{U + D) and xg at Q^ ^ 10 GeV^ 



valence distributions in a previously unmeasured region of x. The early LHC data are not yet 
discriminating. 

Another way to compare PDF predictions is to look at parton-parton luminosities. Fig. [5] 
shows the 2010 and 2011 q — q and g — g luminosity plots. The 2011 comparisons are is 
considerably better agreement. 



4 The value of as{Mz] 



All groups bar GJR use values ~ 0.118 — 0.120 at NLO but there is a definite low(0.113)- 
high(0.117) split at NNLO. MSTW obtain the highest value at both NLO and NNLO and this 
has been atributed to the use of Tevatron jet data in their fits. However, ABM have tried 
inputting these data and find this has only a small effect on their as{Mz) extraction jlO) . 
There is also a 'folk-lore' that DIS data prefer lower values of Ug however both MSTW [13] and 
NNPDF [m have performed DIS only fits in which they find that only the BCDMS data prefer 
low Qfs values. The HERA data actually prefer quite high values. This year we have heard from 
HERA themselves [15|. HERA have input jet data from HI and ZEUS in addition to the HERA 
inclusive DIS data to obtain the HERAPDF1.6 fit. This fit also extends the HERAPDF1.5 
parametrisation to use 14 free parameters (this fit is called HERAPDFl.Sf)- a term is added 
such that the gluon may become negative at low x, Q^ if required (though it does not do so in 
the kinematuc region where data is fitted) and the low- a; valence shape of the d-quark is freed 
from that of the u-quark. When as is freed a value of as{Mz) ~ 0.1202 ± 0.0019 is obtained 
(where the error excludes the scale errors). This fit also has a harder high-x gluon density 
(and corresponding lower low-x gluon density) than the HERAPDF1.5 fit, which brings the 
gluon-gluon luminosity plot for HERAPDF1.6 into closer agreement with MSTW2008. Fig. [6] 
shows the y^ profile of HERAPDF fits with free a^, with and without the jet data, from which 
we can appreciate that as can only be determined with the additon of jet data. This figure 
also shows the PDFs from HERAPDF1.6 free as{Mz) and the q — q and g — g luminosity plots 
for the HERAPDF NLO updates of 2011. 
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Figure 4: Upper plots: the ratio of W^ and W cross sections at the LHC for the PDFs 
considered by the PDF4LHC group. Lower plots: the ratio of Z to W^ + W^ cross sections at 
the LHC for the PDFs considered by the PDF4LHC group. Left-hand side April 2010, right- 
hand side April 2011. Plots from G. Watt http://projects.liepforge.Org/mstwpdf/pdf41hc/2010/l 
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Figure 5: Left hand side: the q ~ q luminosity in ratio to that of MSTW2008 for 
various PDFs. Right hand side: the same for the g — g luminosities. Upper row, 
PDFs 2010: bottom row, updates for CTIO and NNPDF2.1 2011. Plots from G.Watt 
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Figure 6: Top left: x^ scan vs as{Mz) for the HERAPDFl.Sf fit without jet data and for 
the HERAPDF1.6 fit which uses HERA jet data. Top riglit the PDFs resufting from the 
HERAPDF1.6 fit witli free as{Mz)- Bottom row the q — q (left) and g — g (right) luminosities 
in ratio to MSTW2008 for the 2011 HERAPDF updates and for HERAPDFl.O. 

5 NNLO PDFs and predictions for the Higgs 

The Higgs cross section is very sensitive to the value of as{Mz) and the gluon-gluon luminosity. 
For Higgs predictions it is also necessary to consider NNLO calculations. Until this year there 
were only a few NNLO PDFS [TB]. The PDF4LHC recommendation concentrated on the 
MSTW2008 PDF beacuse it is a global fit. However the predictions of the JR, ABKM and 
HERAPDFl.O NNLO PDFs all lie to the low side of MSTW2008 and this was used [l7] to 
cast doubt on the Higgs exclusion limits from the Tcvatron which were based on MSTW. Two 
further issues have arisen related to this. Firstly, Accardi et al [18] have reconsidered deuterium 
corrections for the fixed target data. They find much larger corrections than have usually been 
accounted for, and this results in greater uncertainty in the high-x d— quark which feeds into 
the high-x gluon PDF when it is determined from Tevtaron jet data where the d — g process 
provides a substantial part of the cross section. Secondly, ABM [19] have examined the use 
of NMC data in the global fits. The NMC data on F2 have often been used for the PDF fits 
rather than the cross section data. However, the extraction of F2 relies on assumptions on the 
value of Fl which may not be consistent with modern QCD caluclations. ABM find that using 
F2 rather than the cross section raises their extracted values of as erroneously. However both 
MSTW[20] and NNPDF[2T] have repeated this analysis and do not agree with these conclusions 
(at the time of writing this is still not resolved) . 

Since the gluon PDF is so important for the Higgs predictions another issue which has been 
raised is the goodness of fit of the ABKM and HERAPDFs to the Tcvatron jet data. Watt 
and Thorne [20] obtain poor x^ for these data when fitting to the HERAPDFl.O, 1.5 and the 
ABKM09 PDF predictions. However these fits only compare to the central predictions of the 
HERAPDF and ABKM PDFs. A more valid comparison would account for the PDF error 
bands. HERAPDF have input the Tcvatron jet data to their fit and they obtain much better 
X^ {x^/ndp = 1.48 for CDF and 1.35 for DO jets). Significantly the resulting PDFs do not lie 
outside the HERADF1.5 error bands (although they do imply a harder high-x gluon on the 
upper edge of the error band). ABKM have also made their own fits obtaining y^ /ndp = 0.94 



6 



HI and ZEUS HERA I+II PDF Fit 

tf = 10 GeV^ 



HI and ZEUS HERA l+H PDF Fit 





<0.9S 
U 0.9 



HERAPDF1.5NNLO(prei.) 

HERAPDFl.(INNLO(pre[.) 



■ HERAPDFI.5 NNL<) (prel.l 
HERAPDFI.O NNL<) (prel.) 



Figure 7: Top row:HERAPDF1.5NNLO PDFs compared to HERAPDFl.ONNLO PDFS on log 
and linear x scales. Bottom row: q — q and g — g luminosity plots for HERAPDF1.5 and 1.0 
NNLO in ratio to MSTW2008NNLO. 

for DO di-jet data. 

In 2011 many more NNLO PDFs are becoming available. NNPDF presented a preliminary 
NNLO analysis called NNPDF2.5 (which will be called NNPDF2.1 NNLO) and HERAPDF 
presented a new NNLO extraction HERAPDF1.5NNL0 [22], with full accounting for experi- 
mental, model and parametrisation uncertainties, using the extended form of their parametri- 
sation. An NNLO PDF set from the CT group is also expected soon. Scale differences in 
NNLO heavy quark calculations arc significantly reduced such that NNLO predictions should 
be in better agreement regardless of the choices made for these schemes. Fig. [7] compares 
the HERAPDF1.5NNL0 with the preliminary HERAPDFl.O NNLO which was issued only as 
central predictions for two values of a^. The HERAPDF1.5 NNLO fit has a harder high-x 
gluon. The figure also shows the q — q and g — g HERAPDF NNLO luminosities in ratio to 
MSTW2008, illustrating much closer agreement with MSTW2008 for HERAPDF 1.5 than for 
HERAPDFl.O. This, added to the fact that HERAPDF now recommend that a central value of 
as{Mz) = 0.1176 be used at NNLO, brings the Higgs predictions from HERAPDF into much 
closer agreement with those of MSTW2008. 



6 Comparison to early LHC data 

Fig. Elshow comparisons of the HERAPDF1.5 NLO predictions to the early LHC data on the 
W asymmetry from CMS and ATLAS. Predictions from other PDFs show a similar level of 
agreement. This figure also show comparisons of various PDFs to the ATLAS inclusive jet 
data. The W data have also been used for very preliminary evaluations of their impact on the 
PDF uncertainties by both the HERAPDF and the NNPDF groups. NNPDF [23] have used 
a rcwcighting technique to asses the impact of these data, finding improvements of the order 
of ~ 40% in the u and d-quark. HERAPDF have fitted the asymmetry data from CMS in 
addition to the inclusive HERA data and the CDF ZO and W-asymmctry data. They find a 
modest improvement for the valence quark densities at low- a; - the region that Tevatron data 
do not reach. 

Plots illustrating these improvements and further comparisons of the PDFs to each other 
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Figure 8: Top: HERAPDF1.5 predictions for LHC Ty-lepton asymmetry data from CMS and 
ATLAS data. Bottom ATLAS jet data in the central and forward regions in ratio to the 
predictions of CTEQ6.6 and compared to other PDF predictions. 



and to Tevatron and early LHC data can be found at 

https://wiki.bnl.gOv/conferences/images/l/la/Plenary.Cooper-Sarkar.0415.talk.pdf 
See also the review of Watt [M| for recent work on benchmarking the PDFs. 
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